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SOME EFFECTS OF FUSELAGE FLEXIBILITY ON 
LONGITUDINAL STABILITY AND CONTROL 
By Bernard B. Klawans and ^rold I. Johnson 


SUMMARY 


Some effects of fuselage flexihility on longitudinal stahility 
and control of a large aii^plane are determined hy application of a 
semirigid analytical technique. The dynamic characteristics are examined 
by studying the period and damping of the short-period longitudinal modes 
of oscillation of the fuselage and airplane. The static characteristics 
are examined by studying the longitudinal stability margins and the 
elevator -control deflections reqxiired for balance in straight or steady 
maneuvering flight . 

The consideration of fuselage flexibility in the airplane studied 
appears to introduce no serious problems insofar as dynamic longitudinal 
stability in the subcritical speed range is concerned. It is indicated 
that, if desired, future designs may Incorporate somewhat more flexible 
fuselages than that studied and still have dynamic longitudinal charac- 
teristics approximately eqtial to those predicted by quasi-static theory 
and roughly eqxiivalent to those predicted by rigid-body theory. 

The changes in static characteristics caused by increased fuselage 
flexibility appear, in general, as decreased strai^t-flight and increased 
maneuvering-flight stability margins and increased elevator-control deflec- 
tions required for balance. The need for a means of longitudinal balance 
in steady maneuvering flight other than elevator dr horizontal -tail 
deflection is indicated if flight is desired at low altitude and high 
speed with airplanes that have f\iselage natural frequencies much below 
those representative of current design practice. 


INTRODUCTION 


Several modern hi^-speed airplanes incorporate structural com- 
ponents that are relatively more flexible than those previously 11801. 
There has been concern that the increased flexibility might appreciably 
modify the dynamic stability characteristics as predicted by rigid- 
airplane theory. (See, for example, ref. 1 .) Particular concern was 
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felt for the possibility of interaction between structural and stability 
vibratory modes because the natural frequencies of the major structural 
components are approaching the natural frequencies of the short-period 
stability mode. 

The problem of the effect of wing flexibility on dynamic longitu- 
dinal stability of large airplanes has been treated analytically in 
the subcritical speed range by a simplified semirigid method in refer- 
ence 2 . In the present paper, the effects of fuselage flexibility are 
studied by the same method for the same class of airplanes. The major 
aspects of this study include 

(1) lynamic stability as described by period and damping of the 
fuselage and airplane longitudinal modes of. oscillation 

(2) Static stability as described by longitudinal stability margins 

(3) Static longitudinal control as described by elevator-control 
deflections required for steady strai^t or maneuvering flight 

Implications of the resiilts with regard to the design of future air- 
planes of the class considered are disciassed. 


smeoLS 


Cp 

% 

Cm 

Cmo 

c 

D 

E 

El 


fviselage force coefficient, positive iipward 
lift coefficient 

straight-flight balance lift coefficient 

pitching-moment coefficient, positive nose upward 
airplane pitching-moment coefficient at zero lift 


mean aerodynamic chord, 
differential operator. 


ft 

a/a(^) 


energy, ft-lb 

bending stiffness, Ib-ft^ 
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G 

g 

H 

h 

h(x) 


K 


• / \ o 

recurring group of terms, apfccDf j which for a given 

•A 

airplane and altitude varies as the square of o>f 

' acceleration due to gravity, 52.2 ft/sec^ 

nondimensional "bending deflection at tail, ht/c 

deflection due to bending, positive downward, ft 

shape of fuselage deflection curve which for this paper is 
assumed to be (x/x-fc)^ 

moment of inertia about Y-axis, ft-lb-sec^ 

longitudinal stability margin, fraction of c, positive 
when stable 

radius of gyration about Y-axis, /SC, ft 


M 






lift, positive upward 
tail length, ft 

mass, slugs 
mh(x)dx, slugs 


Mg = 


'^n 


mh(x)x dx, slug-ft 




M7 


= m[h(x)]‘ 

Jxt 


‘dx, slugs 


Mq free -stream Mach number 

m mass per running foot, slugs /ft 

Q generalized force or moment 
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d.0 

generalized coordinate; or pitching velocity, — 
free-stream dynamic pressiare, Ib/sq. ft 

wing area, sq. ft 

time, sec 

true airspeed, fps 

stability coordinate axes (defined in fig. 1) 
earth coordinate axes 

distances along stability coordinate axes, ft 
distances along earth coordinate axes, ft 

nondimensional distance between 0.25c and center of gravity, 
c lengths, positive when center of gravity is rearward 
of 0.25c 

angle of attack, positive nose upward, radians 

elevator deflection, positive when trailing edge is 
downward, radians 

downwash angle, radians 

angle of pitch, positive nose vnpward, radians 
aii^lane relative -density coefficient, Mj^/pS^c 

density, slugs/cu ft 

elevator effectiveness factor, Ct /Ct 

Se/ 

fuselage natural frequency, cps 

effective fuselage frequency, 2 K(Df ^ - 
radians / sec 
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Subscripts ; 

A entire aiip)lane 

eg center of gravity 

e effective 

F flexible 

f fuselage 

ff fuselage forward of center of gravity 

fr fuselage rear of center of gravity 

k kinetic 

n nose of airplane 

p potential 

R rigid 

SF straight flight, flexible fuselage 

SR straight flight, rigid fuselage 

TF txnrning flight, flexible fuselage 

TR turning fli^t, rigid fuselage 

t tail of airplane 

typ typical current design practice 

w wing 

X X. stability axis 

Dots are used to indicate differentiation with respect to time; 

for example, d = — . 

dt 


The subscripts a. Da, q, H, DH, and 6© indicate differentia- 


tion; for example, with respect to the lift parameter, = 


SC-r 

^^Da ~ 

2V 


SCl 

2V 


Cr„ = 


SCt 


% Be ^ 


SCt 

— , and 0 
2V 


L8. 


^5e‘ 
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All coefficients, except those designated lay sxiitahle subscripts, 
refer to the entire airplane. 


AUALYSIS 

Basic Considerations 


Some effects of fiiselage flexibility on the longitudinal character- 
istics of a large airplane are examined by means of the Lagrangian 
formulation of the differential equations of motion. Three degrees of 
freedom are allowed: namely, angle of pitch 9, vertical displacement z, 

and fuselage bending h. (See fig- 1.) 

Assumptions . - Unless otherwise stated, the following assumptions are 
made throi:ighout the paper: 

(1) The airplane forward velocity is constant. 

(2) All drag forces are negligible. 

(3) All lift forces on the fuselage are grouped in the wing 
coefficients . 

(it) Wing and tail mass and lift forces are concentrated at their 
respective quarter-chord fuselage stations. 

( 5 ) Wing and tail are considered rigid. 

(6) The fuselage is perfectly elastic. 

( 7 ) The fuselage deflects parabolically regardless of applied load 
distribution - 

(8) Small classical perturbations are allowed. 

Lagrangian formulation .- 
tion made is 

A. ^ 

dt 8q 


he Lagrangian eq^uation within the assimp- 


8c[ 8q 


( 1 ) 


This equation requires only that the energies (both kinetic and potential) 
and the generalized forces of the system under study be known. 
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Energies . - Figure 1 shows that for small perturbation the vertical 
displacement at any fuselage station is 


z'f = z'cg - x0 + h(x)h^ 


( 2 ) 


where 


X X' Z Z' cos 9 as 1 


Differentiation with respect to time yields the vertical velocity 


z'f = z'cg - x0 + h(x)h^. 


( 5 ) 


Thus f the kinetic energy of • the entire airplane is 


Ek = 



|n(z'f)2dx 




and the potential energy is 



mgz’f dx + 



dx2 I 


( 5 ) 


where the first integral represents the stored energy due to vertical 
position and the second integral^ which may be calculated according to 
the method of appendix A, represents the stored energy due to the 
elastic deformation. 
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Generalized forces . - The generalized forces are "the aerodynamic 
forces on the wing and tail which, for small perturbations, may he 
assumed to act vertically. The forces work when they act through a 
distance. Therefore, for a displacement involving all three .degrees of 
freedom. 


A¥ork = - |AZcg - x^0 + [h(x)] - (Az^g - x^AG + Ah^) I»t (6) 


Thus, the generalized forces for each degree of freedom are 


„ _ AWork 


= -Lw - Lt 


( 7 ) 


% 


Ai-7ork 

A0 


Vv + ^t^t 


(8) 



AWork 

Aht 


-&(x)],,Iv - Lt 


( 9 ) 


where the lift on the wing and tail, respectively, are 


Iv - *lo^ “ °w‘lo^ 
^■T c% 

and 

Lt = 


Expressions for and <x^, the angles of attack of the wing and tail 

in terms of the three degrees of freedom allowed, are derived in appendix B. 


1 
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Equations of Motion 

Four sets of differential eqviations of motion were used in describing 
the effects of fuselage flexibility on longitudinal stability and control; 
the semirigid, rigid, quasi-static, and static -balance equations. 

Semirigid forms . - Performing the operation of equation (l) on equa- 
tions (2) to (11) and then converting to nondimensional notation (see 
ref. 2) gives the semirigid forms of the equations of motion for the three 
degrees of freedom allowed: namely, 

vertical force 


2pD(ot> — 0) + 2p + oCt + Dqi^Ct + DSICt + 


M 


La -^2 Loa " 




( 12 ) 


pitching moment 


2|i ^ D^e - 2p ^ D^H - otCjj, - Dcxfejj. - 

c M^c “a 2 “Da 2 ^ 


Dek;„ 


HC^nn - = 0 


(15) 


and fuselage bending force 


2|i ^ D(a - 0 ) - 2ii ^ D^e + 2p D^H + oCp + Do^p^ + 

Ma M^c o 


Ma 


■a 2 -^Da 


DelOp^ - H(Cp^ ^ - =1^ I ■ 0 


(14) 
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The generalized masses (e.g., Mi) in equations (12) to (l4) are 
recirrring constants and are derived and defined in appendix C . The 
eq^uations used for « evaluating the aerodynamic coefficients are given in 
appendix D. 

Rigid forms.- A special case of the equations of motion^ the famil- 
iar rigid-body case^ is found when the fuselage bending is zero. Equa- 
tions ( 12 ) to (l4) then reduce to 

vertical force 


2^D(a - e) + oCl^ + + De|Cj. + - 


Ct =0 


(15) 


and pitching moment 


D^0 - ~ “ ^eQuig " ~ ^ 


Quasi -static solution.- Another solution to the equations of motion 
of general interest^ called the quasi-static solution^ is found when the 
effects of velocity and acceleration in bending on the forces in the 
system are zero. Equations (12) to (l^) then reduce to 

vertical force 


2|iD(a - 0 ) + ocCt 


+ Dak: 


a 


T + D0-;^t 

4)a 2 






- - 0 


(17) 


pitching moment 

D^e - oCjj^ - Da|CmDa " " ^eCms^ " CmQ = ° 
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and fuselage bending force 


S 'tall c -balance f omi . - Still another form of the semirigid equations 
of motion, which describes static balance or trim of an airplane with a 
flexible fuselage in steady maneuvering flight, is obtained when the rate 
of change of angle of attack and the pitching acceleration as well as the 
velocity and acceleration in bending are all set equal to zero. Equa- 
tions (12) to (l4) then reduce to 

vertical force 


-2pD0 + oC]^ + + HCi^ + 5eCLg^ - Clq = 0 (20) 

pitching moment 


- Cjjjq - 0 - (21) 


and fuselage bending force 


-2nD0 ^ + SaOFs^ " °IO ^ = 0 <22) 


Equations (20) to (22) describe the static balance or trim of an 
airplane with a flexible fuselage in steady straight flight if all the 
terms involving D0 are deleted. 
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Dynamic Longitudinal Stability 

An exponential solution to the equations of motion was assumed. If 
such a solution exists, the d.eterminant of the coefficients must vanish. 
This determinant when written as an equation yields what is commonly 
called the characteristic equation of the system. The roots of the 
operator D contained in the characteristic equation determine the period 
and the time to damp of the various modes of motion in the system mder 
stu(^ if the mode is oscillatory or simply the rate of convergence or 
divergence if the mode is nonoscillatory. 

The foregoing procedure is described in detail in reference 3 
other literature and was used herein to determine the dynamic longitudinal 
stability characteristics of the semirigid, rigid, and quasi-static air- 
plane systems considered. 


Static Longitudinal Stability and Control 

The effects of fuselage flexibility on the static longitudinal sta- 
bility and control are examined by comparing the classic straight-flight 
and maneuvering-flight stability margins and elevator angles necessary 
for balance for rigid airplanes with those for airplanes in which fuse- 
lage bending occurs. 

Stability margins . - The straight-flight stability margin is defined 
herein as the longitudinal distance expressed as a fraction of the length 
of the wing aerodynamic chord between the center of gravity and the point 
at which the lift associated with incremental lift coefficient acts;^ the 
incremental lift coefficient is produced by an infinitesimal change in 
steady forward flight speed along a straight and level flight path with 
controls fixed in the position required for balance at the particular 
initial speed under consideration. VThen expressed analytically, the 
straight-flight stability margin is 


Ks = 


SClq^SV 


(23) 


The static^halance equations (20) to (22) may he differentiated vith 
respect to velocity to obtain and bCj^/bY where the pitching 


^It should be noted that this point is not necessarily a fixed point 
on the airframe in the case of a flexible airplane but^ rather^ is depen- 
dent on center-of -gravity location and stiffness characteristics. 
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velocity D9 may be neglected in level xmaccelerated flight. Substitu- 
tion of these values into equation (25) yields the straight-flight sta- 
bility margin for an airplane with a flexible fuselage as 




Cm G 

y2 

Ct„C-r> - Ct. C-n>_ - C 






a 


7 “= 






a 


— Ct Ctii 

I'a % 



( 24 ) 


which, for infinite stiffness (g = oo)^ reduces to the well-known rigid- 
body straight-flight stability margin 


%R = - 



( 25 ) 


The maneuvering-flight stability margin is herein defined as the 
longitudinal distance expressed as a fraction of the length of the wing 
mean aerodynamic chord betvreen the center of gravity and the point at 
which the lift associated with incremental lift coefficient acts; the 
incremental lift coefficient is produced by an infinitesimal change in 
steady angle of attack at constant forward velocity with controls fixed 
in the position required for either the initial straight or the initial 
steady curvilinear flight path at the flight speed under consideration. 
V/hen expressed analytically, the maneuvering-flight stability margin is 


% = 




(26) 


. The static -balance equations (20) to (22) may be differentiated with 
respect to angle of attack a to determine the solution of the maneuvering- 
flight stability margin Kqi as defined'by equation (26). For an airplane 
with a flexible fuselage, becomes 
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M- 


Kijp - 




a“H 


+ Cjjj Cp_ + c. 




v2 




2^^(CLHCF^ - Cl^% - Cl^ 


( 27 ) 


which, for infinite stiffness (G = <»), reduces to the well-known rigid- 
hody maneuvering-flight stability margin 




1^ 

2n 


(28) 


The stability mar gins presented herein for airplanes with flexible 
fuselages represent real stability rather than apparent stability which 
normally is measured directly in flight tests from the variation of 
elevator angle with speed or acceleration. In this connection it should 
be noted that, for an airplane with a flexible fuselage, some of the 
change in elevator angle needed to change speed or acceleration is 
absorbed by fuselage deflection so that the elevator deflection measured 
■in flight tests cannot, in general, be used as a measure of real stability. 

Elevator -control characteristics in straight and steady maneuvering 
flight.- Thft static -balance equations (20) to (22) may be solved for 
elevator angle necessary for balance in level Ig flight where the pitching 
velocity D9 is 0. Thus, for an airplane with a flexible fuselage. 




CLq 


'SF 


c%2 


(cLhCf^ - Cl^Cf^ - Cl^ 

C-r (c-m - Ct — 
^6eV “a c iv2 


( 29 ) 
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whicli, for infinite stiffness (G = co)^ reduces to 


^esR - 




( 50 ) 


The elevator -angle increment necessary to achieve steady flight 
different from Ig is obtained from the static -balance equations (20) to 
(22) for D0 5^ 0 with all Ig terms neglected. Thus, for an airplane 
with a flexible fuselage. 


A5 






^6er“a 


Ct. 


t\ G 


Ee 


%e K - % im 


a 


c /y2 


( 31 ) 


which, for infinite stiffness (G = oo)^ reduces to 



(32) 


It should be noted that any discussion pertinent to the effects of 
zero-lift pitching moment on static stability or control is based 

on the assumption that the wing and tail are rigid. 
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NUMERICAL CONSTANTS 


Several constants were chosen to identify the study with a typical 
current large airplane: 


Iy = 1.12 X 10^ lh-ft-sec 2 
^A 


- 4.927 per radian 

°w+f 


Ct = 0.802 per radian (based on wing area) t = 0.5 
at 


^ = 0.45 

da 


= 2.72 cps c = 13 f t ■ S^ = 1,428 S(1 ft 


The dimensions and mass distribution shown in appendix C were 
selected such that (l) the moment of inertia of the assumed airplane 
approximated 1.12 x 10^ lb-ft-sec2 over the center-of -gravity range 
investigated (0.25c to 0 . 544 c), and (2) the mass distribution was 
describable analytically. This second criterion reduced the evaluation 
of the mass integrals from a tedious graphical method to one in which 
the values are determined by expressions obtained from simple analyti- 
cal integrations . (See appendix C . ) 


RESULTS AND DISCUSSION 


The effect of fuselage flexibility on the characteristics of a 
bomber airplane is discussed with regard to 

(1) Dynamic stability as described by period and damping of the 
fuselage and airplane longitudinal modes of oscillation 

(2) Static stability as described by longitudinal stability margins 

(3) Static longitudinal control as described by elevator-control 
deflections needed in steady straight or maneuvering flight 


Period and Damping 

Aerodynamic coefficients, mass distributions, and a range of center- 
of -gravity positions were chosen to represent a current high-speed bomber 



NACA W 55^5 


17 


airplane. These values, shown in table I, were substituted into equa- 
tions (12) to (19) to determine period and damping of the airplane and 
fuselage modes of oscillation at altitudes of 8,000 and 50,000 feet. 

(See ref. 5 the method.) The results are summarized in figure 2 
over the entire possible fuselage natural -frequency range at a Mach ■ 
number of 0.7* Figure 3 summarizes the period and damping of the airplane 
mode as a function of dynamic pressure for a center-of -gravity location 
at 0.255 for foiar representative fuselage natural frequencies. The fre- 
quencies selected were (l) oif = cxj or rigid body, (2) ojp = 2.72 cps 

typ 

or typical, current bomber practice, (5) o>f = I.36 cps or one-half the 

frequency of cof , and (it) cdp = 0 or completely flexible body, 
typ ^ 


Rigid-body solution .- The rigid-body characteristics as obtained from 
equations (15) and (16) may be identified in the figures as those at 
oif = ". In figure 2, the center-of -gravity location appears to have little 

effect on the time to damp but causes large changes in the period of 
oscillation. The increase in period with rearward movement of the center 
of gravity is caused by the decrease in restoring pitching moments asso- 
ciated with the reduced stability margins . When the center of gravity is 
at the static neutral point (0.5^45), the motion consists of two aperiodic 
modes, both of which have positive damping . 

Airplane mode of the semirigid body .- The airplane short-period mode 
of oscillation appears to have reduced damping as the fuselage becomes 
more flexible. This is probably caused by the reduced tail rotary damping 
which becomes zero at = 0. However, the airplane mode is still damped 

at dDjp = 0 because the vertical motion of the wing and fuselage supplies 

damping of the plunging motion involved in the short-period oscillation. 
Becaiise of the assumptions made in the analysis, the wing motion contri- 
butes no rotary damping when the center of gravity is at 0.25c and = 0 

but, as the center of gravity moves rearward, some rotary damping due to 
the wing is probably introduced and may account for the increased damping 
shown by figure 2. 

The decrease in period of the airplane mode with increasing fuselage 
flexibility is believed to be due to an apparent forward movement of the 
effective center of gravity as the fuselage becomes more flexible, ^/hen 
the fuselage becomes completely flexible, the airplane probably tends to 
behave like two "tailless" aiiplanes hinged together flying in tandem. 

Each tailless airplane of the combination woiild have a relatively large 
stability margin and much smaller moment of inertia than the origiml 
airplane. It might be expected that the airplane mode characteristics 
would be determined mostly by the largest and heaviest tailless airplane, 
which would consist of most of the fuselage and the wing of the original 
conventional airplane. 
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Fijselage mode of the semirigid body .- Inclusion of a third degree 
of freedom, that of fiiselage flexibility, introduces a new mode of Oscil- 
lation. TMs mode represents the motion the flexible f\iselage would 
assinne. Although not shown in figure 2, the period of the rigid fuse- 
lage = w) is zero because the structural restoring forces are infinite. 

As the fuselage natural frequency decreases, the structural restoring 

forces decrease and disappear at = 0. Existence of a fuselage oscil- 
lation at = 0 (see fig. 2) therefore indicates that a restoring force 
still exists that may be attributed to the air forces present. 

Also, although not shown in figure 2, the time for the fuselage 
oscillation to damp to 0.1 amplitude for oif = “ would be Infinite because 

structural damping is not considered, and the effect of air damping wovild 
be negligible in comparison with the infinite struct\aral restoring forces. 
As ccf decreases, air damping becomes important and causes a gradual 

increase in damping as shown in figure 2. At oif = 0, the damping as well 

as the period of the fuselage oscillation is determined entirely by the 
air forces. 

Quasi-static solution .- Characteristics based on eqiiations (17O to 
(19), wherein the effects of velocity and acceleration due to bending on 
the forces of the system are neglected, are much more readily obtained 
than the semirigid solutions since the characteristic equation is reduced 
from the fourth to the second order. Agreement between the two solutions 
appears to be excellent (see fig. 2) except for exceedingly low fuselage 
natural frequencies for which the quasi-static method predicts'' too much 
damping. 

Effect of altitude .- The rigid-body characteristics indicate that 
both the period and the time to damp to 0.1 amplitude increase with alti- 
tude at constant Mach number (fig. 2). This increase is due to the 
increased ratio between the airplane mass forces and the aerodynamic 
restoring and damping forces which is brought about by a reduction in 
dynamic pressure at the higher altitudes. 

The period of oscillation is essentially independent of altitude at 
any given dynamic pressure (fig. 3 )« However, the time to damp to 0.1 
amplitude at a given dynamic pressure increases with altitude because the 
increase in true speed leads to smaller tail and wing angles of attack 
because of pitching velocity and, hence, smaller rotary damping moments. 

The inclusion of flexibility does not seem to alter these basic trends. 

Significance of flexibility .- It is apparent from either figure 2 or 
figure 3 that fxiselage flexibilities associated with natural frequencies 
above the order of 1 cps do not markedly change the period of damping from 



NACA TN 55^5 


19 


that ejq)ected with a rigid, fuselage . Thus, within the framework of the 
present analysis , designers may possibly incorporate fuselages with 
appreciably lower natural frequencies than that of the typical fuselage 
chosen (2.72 cps) without great concern as to the adverse effects on the 
dynamic longitudinal stability characteristics of the proposed airplane. 

However, when the period of the short-period oscillation approaches 
that of the fiiselage mode which occurs in flight as oif — >0, there is some 

evidence of adverse interaction between structural and stability modes. 
This adverse interaction is indicated by the rapid increase in time to 
damp to 0.1 amplitude as — ^0 shown by the results obtained with the 

semirigid method. It should be noted that the quasi-static method, which 
does not admit the existence of a fuselage oscillation, does not show a 
corresponding rapid increase in time to damp to 0.1 amplitude as oif — ^0. 

The aerodynamic restoring forces which act on an airplane, as 
pointed out previously, are Indicated to be sufficient to increase the 
frequency of the f^lselage mode in flight so that the latter frequency 
is always higher than either the fuselage ground structural frequency 
or the airplane -stability natural frequency. Therefore, the probability 
of instability associated with proximity of the -frequencies of the fuse- 
lage and airplane modes is minimizedj in the present case, indeed, there 
was no indication of instability due to interaction between the two 
modes . 


Longitudinal Stability Margins 


Straight- and maneuvering-flight stability margins were calculated 
according to equations (2h) and (27) and are presented in figmres 4 and 
5, respectively. The results are given as the ratio of the stability 
margin of an airplane flying at a Mach number of 0.7 at 8,000 and 
50,000 feet with a flexible fuselage to the stability margin of this 
airplane operating under identical conditions but with a rigid body. 


Straight flight .- The straight-flight stability margin of an air- 
plane that has zero pitching moment at zero lift = Oj is reduced 

as fuselage flexibility is increased and, in fact, ^becomes zero for 
an airplane with an infinitely flexible fuselage. (See fig. 4.) Center- 
of -gravity location over the range investigated has little if any effect 
on the stability -margin ratio at a given fuselage frequency. 


However, the situation is markedly changed when there is a finite 
pitching moment about the wing quarter chord at zero lift ^ 0^. 
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For a value of C~ given "by 


Cfflo = - -z — Clq 


the stability margin for a flexible fuselage remains constant over the 
entire range of flexibility from an infinitely rigid to an infinitely 
flexible body. Any value of above this shows an increased stabil- 

ity margin with increasing flexibility, and, conversely, any value below 
this shows a decrease in stability margin with increased flexibility. 

The center-of -gravity location has an effect on the degrees of sta- 
bility increase or decrease when Cju^ ^ 0. When the center of gravity is 

at the rigid-airplane neutral point, the ratio of stability margins is, 
of course, infinite (either positively or negatively depending on the 
sign of Cdiq) with respect to the rigid-body margin. At more forward 

center-of -gravity locations, the ratio of stability margins approaches 
the values for Cjuq = 0. Therefore, it is seen that the importance of 

a given value of Cjjjq is inversely proportional to the initial rigid- 

airplane stability margin. 

Increasing altitude at constant Mach number reduces the effects of 
flexibility on the static margins for any airplane configuration. This 
effect is believed to be solely due to the relative reduction in air 
forces compared with the structural forces because of the reduced dy n a mi c 
pressure at any constant Mach niomber that occurs with increasing altitude. 

Maneuvering flight .- As the fuselage becomes flexible , the 

maneuvering-flight stability margin shows a marked increase over the 
rigid-body values. Although not readily seen from figure 5^ "tlie increased 
margins amount to about 20 percent chord over the center-of -gravity range 
investigated at = 0. The increase in maneuvering-flight stability 

margin ratio with increased altitude probably stems from the reduction in 
rigid-airplane stability margin which normally occurs with increasing 
altitude and is associated with reduced damping in pitch. Because the 
turning-flight stability margins involve no change in speed, the zero- 
lift pitching-moment coefficient has no effect on these margins even 

when flexibility is considered. 
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Significance of flexibility »- The conclusion often reached is that 
increasing the flexibility of a fuselage will always result in decreased 
stability margins because of the tendency of the tail to stream with the 
relative wind and thereby produce less lift than the same tail on a 
rigid fuselage. Such an effect is apparently predominant when straight- 
flight stability margins are considered. However, although such a 
streaming tendency exists in steady maneuvering fli^t, the present analy 
sis indicates that the effect of normal acceleration on the distributed 
fuselage mass for airplanes with conventional tails causes both ends of 
the fuselage to deflect downward siafficiently such that the tail- of the 
flexible airplane produces more lift than that of the rigid airplane. 

The resulting increase in tail lift draws the centroid of the total air- 
plane lift rearward which is synonomous with increased stability margin. 
Thus, although the straight -flight stability margin goes to zero 
^for Cjj^^ = 0^ as the fuselage flexibility is increased, the maneuvering- 

f light stability margin is increased. 

For flexible airplanes, the problem of straight-flight stability 
margin may be minimized through carefiil design practice. Advantage can 
be taken of the situation when 



If the wing sections are selected such that this eq_uation is obeyed at 
the design speed and altitude condition, flexibility will have no effect 
and the airplane will have the same stability it would have if it had a 
rigid fuselage . As the airplane operates off the design-speed conditions 
the effects of fuselage flexibility become important. 


Except for large negative values of which in most cases can be 

avoided in airplane design, the effects of the inclusion of fuselage- 


flexibility effects on current airplane design 




typ 


= 2.72 cps) with 


regard to stability margins are either negligible or beneficial in both 
straight and maneuvering flight. Indeed, from figures 4 and 5^ designers 
may expect to use fuselage natural frequencies appreciably less than is 
the current practice without expecting very large differences. 


Elevator -Control Motion 


The elevator angles necessary for balance in either straight or 
maneuvering flight were calculated according to equations ( 29 ) and ( 5 I) 
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and are presented in figures 6 and J , respectively. In the case of 
straight flight, the elevator angles necessary for balance are given as 
absolute values. In the case of maneuvering flight, elevator-angle 
increments necessary for balance at any given normal acceleration dif- 
ferent from 1 g are given as the ratio of the elevator -angle increment 
req.xiired for an airplane with a flexible fuselage to the elevator -angle 
increment req,\aired for an airplane with a rigid body. In every case, 
the airplane is flying at a Mach niimber of 0-7 at either 8,000- or 
50 , 000 -foot altitudes. 

Rigid-body control .- The elevator angles at ojj. = « in figiore 6 
represent those expected from a rigid airplane of the assima.ed geometry. 
The static stick-fixed neutral point occurs when 6g = 0 for Cmo “ 

which is at the center-of -gravity location of 0.5^^c. In this condition, 
a positive lift necessary for balance is generated on the entire tail 
when the rigid airplane is at an angle of attack. As the center of 
gravity moves forward, less tail lift is req_iiired for balance. This is 
achieved by increasing up-elevator (negative) angle as may be seen in 
figure 6. IThen the center of gravity reaches the wing q.i;iarter-chord 
position, the up-elevator angle required is just sufficient to eliminate 
completely all lift on the horizontal tail for the configuration which 
has no zero-lift pitching moments or distributed fuselage lift. Of 
course, positive Cjjj^ causes a down-elevator angle necessary for balance 

and negative Cjjj^ causes an up-elevator angle necessary for balance with 

respect to that required when Cjjjq = 0.. 


Straight-flight control of flexible fuselage . - The elevator angle 
necessary for balance in straight flight for Cjjj^ = 0 becomes increasingly 

negative (trailing edge up) as fuselage flexibility is increased and, in 
fact, becomes infinite for an airplane with an infinitely flexible fuse- 
lage. (See fig. 6.) However, when zero-lift pitching moments are con- 
sidered, it can be shown that fuselage flexibility has no effect on the 
elevator angle necessary for balance when 





+ Cm c. 






Ct—C-r' - C]^ Cp 




■a -^H 


More positive values of Cjjjq reqxiire am increasing positive elevator 

angle until as aj- — ^0 the elevator angle necessary for balance becomes 
infinite. Conversely, for more negative values of C^^, the elevator 

angle necessary for balance becomes infinite in the opposite direction. 
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Increasing altitude requires increased up-elevator angles necessary 
for balance. This is probably due primarily to the reduced dynamic pres- 
sure which requires that the airplane be flown at higher lift coefficients 
at high altitudes . 

Maneuvering -flight control of flexible fviselage .- The maneuvering- 
f light eleva tor -angle -increment ratio necessary for balance for an 
airplane with a flexible fuselage increases as fuselage flexibility is 
increased and becomes infinite for an airplane with an infinitely flexible 
fuselage. (See fig. 7.) 

Significance of flexibility . - The deflection curve which a flexible 
fuselage assumes depends on the distribution of applied aerodynamic loads, 
the mass distribution, and the stiffness distribution. In no rmal flight, 
that is, with the center of gravity at or near O.256, the air load 
required on the horizontal tail for balance in straight flight is zero if 
Cj^ = 0. Und.er these conditions, the fuselage may be tho\ight of as a long 

flexible beam having distributed mass which is supported by a single con- 
centrated lift force supplied somewhere near the middle of the beam. The 
result is that both ends of the fuselage will droop, and it is apparent 
that increasing up elevator will be required to preserve balance as the 
fuselage increases in flexibility. When the center of gravity is moved 
rearward from 0.25c, up tail loads are needed for balance. If the up 
tail load for balance exceeds the weight of the tail assembly (which 
occurs for the airplane studied when the center of gravity is about 0.52c), 
the fuselage deflection curve begins to have an inflection point at the 
rear end which moves fonrard with further rearward movement of the center 
of gravity. Therefore, the assumption of a simple deflection curve such 
as the parabola becomes less satisfactory as the center of gravity 
approaches the neutral point. However, it is apparent from the foregoing 
discussion that the assuirption of a parabolic deflection curve is at least 
qualitatively realistic and should allow some idea to be gained in regard 
to the onset of serio\is trim problems . 

The increases in elevator angle compared with those required for a 
rigid body = 00^ are not marked for the design studied = 2 .72 cps 

for normal center-of -gravity locations (approximately 0.25c) for the air- 
plane in level flight at both 8,000 and 50,000 feet and also for the air- 
plane at 50,000 feet in maneuvering flight. In fdct, from figures 6 
and 7 i'fc appears that designers can reduce appreciably current fviselage 
structural frequencies before any serious elevator-control balance prob- 
lems will be expected to occur. However, the elevator-control balance 
problem in low-altitude high-speed flight (Mq = O.7O) appears in figure 7 
to be becoming serious for current airplanes in maneuvering flight and 
will probably become intolerable for more flexible designs. If high-speed 
flight becomes necessary at low altitudes for airplanes having fviselage 
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natural freq.uencies much, helOT those of present-day practice, some means 
other than elevator or stabilizer control may become necesseiry to balance 
the airplane longitudinally in maneuvering flight. 

It is interesting to note from figinre 6 or figure 7 that the eleva- 
tor is, in general, corapletely incapable of balancing the airplane with a 
completely flexible fuselage in either straight or txurning flight. Thus, 
any discussions of static stability for fuselage stiffness very close to 
zero are primarily of academic interest. 


■ CONCLUDING REMAEKS 


The consideration of fuselage flexibility in a current large high- 
speed airplane design appears to introduce no serious problems insofar 
as dynamic longitudinal stability in the subcritical speed range is con- 
cerned. It is indicated that, if desired, future- designs may incorporate 
fuselages with natural frequencies appreciably less than that studied and 
still have dynamic longitudinal characteristics approximately equal to 
those predicted by quasi-static theoiy and roughly equivalent to those 
predicted by rigid-body theory. 

The aerodynamic restoring forces which act on the airplane are 
sufficient to increase the frequency of the fuselage oscillatory mode in 
flight so that the latter frequency is always higher than either the 
f\iselage ground structmral frequency of the airplane -stability natural 
frequency. Thus, the occiarrence of a resonant condition is avoided. 

The changes in static characteristics caused by increased fiiselage 
flexibility, in general, appear as decreased straight-flight and increased 
maneuvering-flight stability margins and increased elevator-control 
deflections reqtiired for balance. The need for a means of longitudinal 
balance in steady maneuvering fli^t other than elevator or horizontal- 
tail deflection is indicated if flight is desired at low altitude and 
high speed with airplanes that have fuselage natural frequencies much 
below those representative of current design practice. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 23 , 1955* 
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APPENDIX A 


CALCULATION OF POTEINTIAL STRAIN ENERGY 


Nodal line 



The potential strain energy is derived in terms of the ma viTtiimi 
Kinetic energy of a harmonically vihrating free -free beam (the fuselage). 
Mathematically, this equality can he expressed as maximum potential strain 
energy equal to maximum kinetic energy or 


d^[h(x)]xt.] 


The problem then is to express Zf in known terms. 

The vertical displacement and acceleration of any fuselage station 
may be expressed as 


dx = 




mzf dx 


(Al) 





H = 


Zcg - X0 + h(x)lL(. 


(A2) 



26 


NACA TM 55^5 


and 


Z|. = Zgg - x0 + li(x)h^ 


(A5) 


It is a necessary condition that the external vertical force and pitching 
moments of a free -free beam he zero at all times. Hence, 


and 



mzf dx = 0 = 


z'c^A + %^1 


(Ah) 



mxzf dx = 0 = ^^2 


(A5) 


Solving equations (A4) and (A5) for z^g and 0, respectively, in terms 
of yields 


% » 

“ Mfl ^ 


•• Mo .. 

0 = -^ h+ 

■% ^ 


(a6) 


or, combining equations (AJ) and (a 6) gives 




- — X + h(x) 
% 


(A?) 
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Under the assiio5)tion of sinusoidal motion it can he shown easily that 


= — ( 




(A8) 


and 


H = 




(A9) 


so that equation (A7) hecomes 


= ht 


X + h(x) 

% ly 


(AlO) 


Equation (AlO) then is an expression for in known terms. Sub- 

stitution of equation (AlO) into eqmtion (Al) and simplifying yields the 
potential strain energy expressed as 


1 r"- 


El ^ 


d2 [h(x)h^] 

2 

(3x = ^ Mzh+^ (l. 

M^2 


dx2 

>• ox — 1 j. 

M#5 

lyM^y 


(All) 
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APrEDIDIX B 

DERIVATIONS OF TAIL AND WING ANGLES OF ATTACK 
For a rigid airplane, the wing angle of attack is 

_ 

%" V 


(Bl) 


and the tail angle of attack is 




V 


(B 2 ) 


For the flexible aiiplane, cl ^ and are modified because of 

fuselage bending displacement and. velocity; conseq.uently. 


°wp ~ 


and 




[d[h(x)h^ 


h(x)h^ 

dx 

Jx=x^ 

J 

Jd[h(x)lLt] 

V -i- 


1 dx 

/ * 
x=x.|. 

V 


(B5) 


Jx— 


m 


JX=Xt 


For the assumed parabolic fuselage deflection under load. 
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so that 


<3-[h(x)] ^ _2x_ 
dx x^2 


(b6) 


Substitution of eq.uations (Bl) into (B5) and (B2) into (b 4) and 
nondimensionalizing gives, respectively, 


= a 



2 ^ J_ 
xt xt 


H + 



(B7) 


and 


cc^ = a 


^fa 

da\ 


Da 



^ D0 - 2H -^ + DH 
c 


(b8) 


Equations (B7) and (b 8) are used herein to describe, the angles of 
attack of the wing and tail for the several degrees of freedom allowed. 
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APPENDIX C 

DERIVATION OF GENERALIZED MASSES 

The assumed loading for the airplane under consideration may be 
pictinred as 



c = 13 ft 

^A “ slugs (125^000 lb) 

= 1,275 slugs (41,000 lb) 

Mf + M^. = = 2,609 slugs 


M^- = 621 slugs 
Mf = Mff + = 1,988 slugs 
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and where, for the assumed fuselage geometry and triangular load 
hution in slugs. 


Mff = 1,299 - 87 OX 


The mass per running foot is 


and 


m^^ = 2M^^ 



(for X = Xn to x = 0) 


-x+ + X 
“f r = 2Mfr — 5— 

xt 


(for X = 0 to X = xt) 


for 


Xn = 51-3 + 13x 


and 


x^ = -46.5 + 13x 


Proceeding to the mass integrals gives 


M]_ = 




mh(x)dx 


distri- 


(Cl) 


(C2) 


(C5) 


(C4) 


(C5) 
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Conibiiiiiig eq.uations (Cl) to (C5) yields 




(C8) 


(C9) 


Eq.uations (C7) to (C9) are iised to evaluate the mass iiitegrals over the 
center-of -gravity range of the investigation. 
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APPENDIX D 

EQUATIONS FOR EVALUATING THE AERODYNAMIC PARAMETERS 
The aerodynamic parameters in the general case may be expressed as 


f(n)=Lf(n) = 


f(n)C. 






■%(n) “ % (n)°I-a„ T * T 




For exao5)le, the airplane lift-curve slope the damping in 

pitch |cmq, and the bending-force coefficient due to velocity in 
bending ^CFpg escpressed, respectively^ as 


oCl = cxCt + a(l - 

V da; 

D^Fdh = 
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Parameter 


Values of parameter for 
center-of -gravity locations of - 



0.255 

0.356 

0.455 

0.5445 

Mi/Ma 

0,1575 

0.1441 

0.1510 

0.1580 



1.762 

1.756 

1.710 

1.685 

M2/M^c 

-.01117 

-.00577 

.00576 

.01148 

M5K 

.07056 

.07718 

.08453 

.09180 



5.568 

5.368 

5.368 

5.368 



1.292 

1.256 

1.220 

1.186 



2.871 

2.298 

1.725 

1.185 



.4482 

• 3797 

.3021 

.2197 

^Ldh 

.8021 

.8062 

.8193 

.8416 

Cl. 

°e 

.4010 

.4010 

.4010 

.4010 



-1.579 

- 1.042 

-.5060 

0 

Ic 



- 4.624 

-4.370 

- 4.125 

-3.896. 

in„ ... 

^mq .... 

-10.280 

-9.758 

-9.356 

-9.082 



-1.604 

-1.612 

-1.639 

-1.683 

in 

2 ^“DH 

-3.155 

-3.055 

-2.972 

- 2.888 

Cm. 

°e 

-1.865 

-1.813 

-1.762 

-1.712 

C-p. 

. 4412 . 

.4453 

.4584 

.4807 

¥^jkt 

1.292 

1.256 

1.220 

1.186 

¥^q 

2.871 

2.790 

2.707 

2.623 



.4482 

.4610 

.4741 

.4862 

¥^im 

.8021 

.8062 

.8194 

.8419 

• 

.4010 

.4010 

.4010 

.4010 










Semirigid solution 
Quasi-static solution 
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Figure 2.- Semirigid eind quasi-static solutions for period and damping of short-period longittidinal 
mode and fuselage mode of a large airplane as a function of fuselage nattiral frequency for 
various airplane center -of -gravity locations. Ifcich number, 0.7; rigid-body values at = oo. 












Figure 4.- Variation of static-stability-margin ratio with fuselage natxiral frequency for a large 
airplane in straight flight at several airplane center-of -gravity locations. Mach nuniber, 0.7 




Figure 5.- Variation of maneuvering-flight statllity-margin ratio with fuselage natural frequency 
for a large airplane at several airplane center-of -gravity locations. Mach number, 0.7. 















